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ABSTRACT 

We construct orbit-based axisymmetric dynamical models for the globular cluster M15 which fit groundbased 
line-of-sight velocities and Hubble Space Telescope line-of-sight velocities and proper motions. This allows us 
to constrain the variation of the mass-to-light ratio M/L as a function of radius in the cluster, and to measure the 
distance and inclination of the cluster. We obtain a best-fitting inclination of 60° ± 15°, a dynamical distance 
of 10.3 ± 0.4 kpc and an M/L profile with a central peak. The infeiTed mass in the central 0.05 parsec is 3400 
Mq, implying a central density of at least 7.4 x 10® Mq pc“^. We cannot distinguish the nature of the central 
mass concentration. It could be an IMBH or it could be large number of compact objects, or it could be a 
combination. The central 4 arcsec of M15 appears to contain a rapidly spinning core, and we speculate on its 
origin. 

Subject headings: distance scale — globular clusters; individual (M15) — stellar dynamics — stars; kinematics 
— black hole physics 


1. INTRODUCTION 

M15 is a well-studied globular cluster. It has a very 
steep central luminosity profile, and may be in the post-core- 
collapse stage (e.g., Phinney 1993; Trager, King & Djorgov- 
ski 1995). Measurements of nearly two thousand line-of-sight 
velocities (from the ground and with HST) and proper mo¬ 
tions (with HST) have recently become available (Gebhardt 
et al. 2000, hereafter GOO; McNamara, Harrison & Anderson 
2003, hereafter M03). 

McNamara, Harrison & Baumgardt (2004, hereafter M04) 
restricted themselves to the subset of 237 stars inside 0.'3 of 
the center of M15 for which both Fabry-Perot radial veloci¬ 
ties and HST proper motions were measured, and computed 
the mean dispersions in these measurements. Assuming the 
cluster is an isotropic sphere, and the observed stars are rep¬ 
resentative, the ratio of these dispersions (one in kms“*, the 
other in mas yr“’) provides the distance (Cudworth 1979, Bin- 
ney & Tremaine 1987). M04 find a distance of 9.98 ±0.47 
kpc, which is consistent with the canonical value of 10.4 kpc 
(DuiTell & Hands 1993), but is smaller than, e.g., the recent 
determination of 11.2 kpc by Kraft & Ivans (2003) who used 
a globular cluster metallicity scale, based upon Fe II lines. 

Here we extend the M04 study by using a larger frac¬ 
tion of the line-of-sight velocity and proper motion samples, 
and comparing these with more general dynamical models 
to study the internal structure of the cluster as a function of 
radius. We follow the approach taken by van de Ven et al. 
(2005, herafter V05), who constructed axisymmetric dynami¬ 
cal models for the globular cluster w Centauri and fitted these 
to groundbased proper motions and line-of-sight velocities. 
This technique provides the internal dynamical structure as 
well as the inclination of the cluster, an unbiased and accurate 
dynamical distance, and the M/L-profile. Our aim is to derive 
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similar information for M15. We are particularly interested 
in the M/L profile, as significant mass segregation is believed 
to have occurred in the cluster (Dull et al. 1997). The HST 
proper motions have sufficient spatial resolution to study the 
dynamical structure and mass concentration in the center. 

In SectionlJl we summarize the observational data. In Sec¬ 
tion 0 we consider the influence of measurement errors on 
the data, select the stars to be used for the dynamical model¬ 
ing and also study the possible residual systematic effects in 
the observed mean motions. We construct dynamical models 
in Section|4] derive a distance, and investigate the effect of the 
unknown inclination and of radial M/L variations in the clus¬ 
ter. We discuss the dynamics of the central 0.2 parsec of M15 
in Section|3 and summarize our conclusions in Section|3 

2. OBSERVATIONAL DATA 

We discuss, in turn, the surface brightness distribution, the 
line-of-sight velocities, and the proper motion data for M15. 

2.1. Surface brightness distribution 

The surface brightness distribution of M15 has been studied 
in detail by a number of authors (Lauer et al. 1991; Trager et 
al. 1995; Guhathakurta et al. 1996; Sosin & King 1997). Noy- 
ola & Gebhardt (2005) reanalysed archival WFPC2 images 
with a new technique which measures the integrated light 
to determine the surface-brightness distribution in the dense 
central regions of the cluster. They combined this with the 
groundbased profile from Trager et al. (1995). The surface- 
brightness profile is a power-law in the inner arcsecond, with 
slope -0.62 ±0.06, which agrees with the determinations by 
Guhathakurta et al. (1996) and Sosin & King (1997). 

Because the surface-brightness distribution of M15 is only 
slightly flattened and the individual stars in the cluster are re¬ 
solved, it is difficult to estimate the ellipticity and the posi¬ 
tion angle (PA) of the major axis. Guhathakurta et al. (1996) 
measured an average ellipticity of e = 0.05 ±0.04. Determi¬ 
nations of the PA of the photometric major axis as measured 
from North through East vary from 125° (White & Shawl 
1987) to 45° (from a DSS image) at large radii, and is given 
as 60° ±20° by Guhathakurta et al. (1996). Gebhardt et al. 
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Fig. 1. — Smoothed radial dispersion profiles of the line-of-sight velocities and proper motions for various cuts in the velocity error distribution, a) Line-of-sight 
velocities. Red is an uncertainty cut of 5 km s“* or less, black is for 7 km s“*, and blue is for 14 km s“'. Uncertainty cuts below 5 km s“* show no difference from 
the red lines. Every error bai' represents a radial ring containing 60 stars, b) major (solid lines) and minor (dashed lines) axis proper motions for various cuts in 
the velocity error distribution. Red is an uncertainty cut of 8 km s“' or less, black is for 14 km s“*, and blue is for 21 km s“* (for an assumed distance of 10 kpc). 
Uncertainty cuts below 8 km s“* show no difference from the red curves. Every error bar represents a radial ring containing approximately 60 stars. 


(1997) reported the kinematic PA of maximum rotation at 
198° North through East. Here we a dopt this value for the 
PA of the major axis (see Section im . 

2.2. Line-of-sight velocity sample 

The bulk of the line-of-sight velocities come from the com¬ 
pilation by GOO, who reported measurements for 1773 stars 
brighter than B magnitude 16.5 out to a radius of 17'. Their 
data set includes earlier measurements by Peterson, Seitzer 
& Cudworth (1989), Dubath & Meylan (1994), Gebhardt et 
al. (1994, 1995, 1997) and Drukier et al. (1998), as well as 
82 stars in the inner region measured using adaptive-optics- 
assisted spectroscopy (GOO). In addition, line-of-sight veloc¬ 
ities for 64 stars in the inner 4" were obtained with STIS on¬ 
board HST (van der Marel et al. 2002; Gerssen et al. 2002, 
2003). 

The expected number of non-members in this data set is 
very small. The cluster is very dense, so few interlopers are 
expected from chance superposition of field stars. The sys¬ 
temic line-of-sight velocity of M15 is -107.5 ±0.2 kms“' 
(GOO). With an internal velocity dispersion of about 12 km s“' 
in the center, the cluster stars are well-separated in velocity 
from the foreground galactic disk. The bulk of the data is in¬ 
side 4' and the median error of the line-of-sight velocities is 
3.5 kms“'. This is a significant fraction of the line-of-sight 
velocity dispersion, especially at larger radii. 

2.3. Proper motions 

M03 published proper motions for 1764 stars within 0(3 of 
the center of M15, derived from multi-epoch HST/WFPC2 
imaging. The stars range in brightness between 14.0 and 
18.3 mag. Of the 1764 stars only 703 stars brighter than B- 
magnitude 16.5 where kept and used in their analysis. We 
restrict ourselves to this sample of 703 stars. 

M03 derived the proper motions in the classical way, by 
using a reference frame consisting of cluster stars, and then 
modeling the difference between the positions on the first and 
second epoch images in terms of a zero point difference, a 
scale change, a rotation, a tilt, and second-order distortion 


corrections (e.g., Vasilevskis et al. 1979). The correspond¬ 
ing transformation equations were then solved using a least 
squares routine. As a result, the derived proper motions are 
not absolute, but may contain a residual global rotation. They 
also may be influenced by perspective rotation caused Iw the 
space motion of the cluster. We return to this in Section|3 

Figure 1 of M04 shows the histogram of all proper motions. 
The median error of these measurements is 0.12 mas yr“', 
corresponding to about 6 km s“* at a distance of 10 kpc, with 
some errors as large as 0.50 mas yr“'. 

3. KINEMATICS 

We use a Cartesian coordinate system (x',y',z') with z' 
along the line of sight and x' and y' in the plane of the sky 
aligned with the cluster such that the y'-axis is the photometric 
minor and rotation axis of the cluster. The kinematic measure¬ 
ments then give v^i for the line-of-sight velocities in kms“\ 
and Pf/ and p,y> for the proper motions in mas yr“'. To convert 
tlx' and /Tv' into Vx’ and v,.' in kms“', we used v.i' = 4.74D/iv' 
and Vy' = 4.74D^y/, with D the distance in kpc. 

The measurement errors in the line-of-sight velocities and 
in the proper motions need to be taken into account when ana¬ 
lyzing the kinematics of M15. We do this by means of a max¬ 
imum likelihood method, which corrects for each individual 
velocity error, and provides robust estimates of the mean ve¬ 
locities and velocity dispersions in spatial bins on the sky. The 
method is described in detail in Appendix A of V05. 

3.1. Selection 

Figure^ shows the velocity dispersion along the line of 
sight, as a function of radius, for three different selections 
of stars. The blue curve is for all stars with velocity error 
smaller than 14 kms“', the black curve for those with errors 
smaller than 7 kms“\ and the red curve for the stars with 
errors smaller than 5 km s“'. Although we corrected the dis¬ 
persion for the individual measurement errors, the curves do 
not all overlap, suggesting that the (larger) errors are not esti¬ 
mated very accurately. The curves converge once we exclude 
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Fig. 2.— Smoothed line-of-sight mean velocity (top) and velocity disper¬ 
sion field (bottom) for M15 in km s“*. 

the stars with errors larger than 7 kms“\ so we restrict our¬ 
selves to this sample of 1546 stars. The corresponding pro¬ 
file varies between approximately 12 km s“' in the center to 3 
kms“' at about 10' (cf. Figure 12 in GOO). 

Figure [Qr shows the velocity dispersions of the proper mo¬ 
tions in the x' and y' direction, respectively, as a function of 
radius, for three different selections of stars based on the mea¬ 
surement errors, corresponding to 21, 14 and 8 kms“^ for an 
assumed distance of 10 kpc. The radial range covered is only 
0.'23, and hence corresponds to the inner data points of Fig¬ 
ure ct. The profiles for the different error selections are con¬ 
sistent to within the errors of the dispersions. The dispersion 
profiles in the two orthogonal directions appear to differ and 
they are, in fact, inconsistent at the formal one sigma level. 
The overall proper motion dispersion is independent of the er¬ 
ror selection. It is therefore not evident that a selection based 
on measurement error is justified, and we use all 703 proper 
motions. 

Figures|2]and|3show a smooth representation of the mean 
velocity and dispersion fields of the line-of-sight velocities 
and proper motions, respectively. The fields where adaptively 
smoothed by computing at each stellar position the kinematic 
moment for the nearest 100 neighbors using Gaussian weight¬ 
ing with distance from the central star. The Gaussian used for 
the weighting has the mean distance of the 100 stars as its 
dispersion. The resulting smooth kinematic maps correlate 
the different values at different points, but bring out the main 
features of the observed kinematics of M15. 

The line-of-sight velocity maps show significant structure, 
including overall rotation (GOO; and Section |3below). The 
line-of-sight velocity dispersion shows the radial fall-off il¬ 
lustrated in Figure!^. The contours of constant dispersion 
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Fig. 3.— Smoothed kinematic velocity fields of the proper motions in 
kms“* assuming a distance of 10 kpc. Left: kinematics in the x'-direction 
(parallel the major axis). Right: kinematics in the y'-direction (parallel the 
minor axis). Top: Velocity field. Bottom: Corresponding dispersion fields. 


TABLE 1 

Systemic proper motion of M 15 


fioi cos (5 
(1) 

B-S 

(2) 

Reference 

(3) 

-0.3 ±1.0 

-4.2 ±1.0 

Cudworth & Hanson (1993) 

-1.0± 1.4 

-10.2± 1.4 

Geffert et al. (1993) 

-0.1 ±0.4 

0.2 ±0.3 

Scholz et al. (1996) 

-2.4 ±1.0 

-8.3 ±1.0 

Odenkirchen et al. (1997) 

-0.95 ±0.51 

-5.63 ±0.50 

Dinescu et al. (1999) 


Note. — Five determinations of the systemic motion of M15 on the 
plane of the sky. Cols (1) & (2): Systemic proper motion in a, 5 in units 
of mas yr“*. Col. (3): Reference. 


are slightly elongated. The mean proper motion maps are dif¬ 
ficult to interpret, because of the relatively large measurement 
errors. The dispersions in the proper motions are nearly con¬ 
stant over the small extent of the field (cf Figure^). 

3.2. Perspective rotation 

The observed motions contain a contribution from the per¬ 
spective rotation caused by the space motion of M15. The 
systemic line-of-sight velocity is -107.5 ±0.2 kms“* (Geb- 
hardt et al. 1997), but the component of the space motion in 
the plane of the sky is not well-determined. Table [TJpresents 
a summary of the reported space motions for M15. 

Use of eq. (6) of V05 shows that any of the values in Ta- 
ble[3result in contributions to the observed proper motions of 
at most 0.0025 mas yr“' at the edge of the small area where 
we have measurements. This is well below the measurement 
errors, and we therefore ignore it. The contribution to the ob¬ 
served line-of-sight velocities is ±0.1 kms“* at 5' from the 
center if we use the Cudworth & Hanson (1993) value for 
the space motion, which is small enough that it can be ig¬ 
nored. If the more recent determinations by Odenkirchen et 
al. (1997) and Dinescu et al. (1999) are correct, then there 
would be a contribution of ±0.8 kms“'. This contribution 
is still not significant relative to the measurement errors and 
therefore we do not apply any correction for the perspective 
rotation. However, in this case perspective rotation would be¬ 
come important for studies of the kinematics near the tidal 
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radius (21 f5, Trager et al. 1995). If not corrected for, it would 
result in an apparent leveling-off of the line-of-sight velocity 
dispersion, as reported by Drukier et al. (1998). 

3.3. Residual global rotation 

It is possible to correct for the possible presence of residual 
global rotation in proper motion measurements of a globular 
cluster by using the line-of-sight velocities and the assump¬ 
tion of axisymmetry. VOS applied this method with success to 
uj Centauri. We follow the same approach for M15. 

In an axisymmetric cluster, the following relation is valid 
between the mean motion and the mean line-of-sight 
velocity (v^/) at any point (x',/); 

(vj,) =4.74Dtan/(py), (1) 

where (v^/) is in units of kms“', D is the distance in kpc, i 
is the inclination of the cluster, and is in masyr“' (cf. 
Evans & de Zeeuw 1994). We computed the mean value ) 
in spatial bins on the plane of the sky, and similarly for (v^/). 
We find a formal best-fit value which corresponds to an incli¬ 
nation of 59°±12° at a distance of 10±0.5 kpc. 

If the proper motion measurements still contain an un¬ 
known amount of global rotation fl (constant with radius), 
then the observed and should be replaced by 

Ma' = fix'+y'^, My/ = My' -x'n. (2) 

We stepped through a range of values of and determined the 
value for which the scatter around the relation dD was mini¬ 
mized. We find a weak minimum for fl equal to -0.21 ± 0.22 
mas (yr arcmin)“’. The difference between the two cases is 
modest and changes the formal best-fit value of the inclination 
slightly to 54° ±10°. We conclude that any effect of residual 
global rotation is below the measurement errors, and therefore 
we do not correct for it. 

4. DYNAMICAL MODELS 

We construct axisymmetric dynamical models of Ml5 by 
means of Schwarzschild’s (1979) orbit superposition method, 
as implemented by Verolme et al. (2002). The inclusion of 
proper motion data is described in detail by V05, together 
with extensive tests designed to establish the accuracy with 
which the distance and internal structure can be recovered. 
These dynamical models are collisionless, which is not 
necessarily a valid assumption for a dense globular cluster 
such as M15. _ 

We start by constructing a luminosity model (Section f4.1> . 
and compute constraints fro m the observed kinematics binned 
into polar apertures (Section f4.2> . In each aperture, we com¬ 
pare the mean velocities and velocity dispersions to the pre¬ 
dictions of the dynamical model while varying the parameters 
to find a best-fit model. The model parameters are the incli¬ 
nation i, the distance D, the mass-to-light ratio (M/L) values 
in different radial bins and a central dark mass Mjark (Sec- 
tion l4.3> . In Section R"^ we o btain the bes t-fit models and we 
discuss the results in Sections f4.5ll4.6l and l4.7l 

4.1. Luminosity model 

We use the Noyola & Gebhardt (2005) V-band one- 
dim ensio nal surface-brightness profile discussed in Sec¬ 
tion |0] as a basis for our mass model of M15. The profile 
extends to 15'. We parameterized the profile with a multi- 
Gaussian expansion method (Monnet, Bacon & Emsellem 
1992; Emsellem, Monnet & Bacon 1994) by means of the 



Fig. 4.— Top: Radial surface-brightness profile of M15 from Noyola & 
Gebhardt (2005) and the corresponding MGE best-fit model consisting of 
14 Gaussians as a function of radius. Each individual Gaussian is shown as 
well. Left axis shows Lq per pc^. Bottom: The residual between the surface 
brighness and the MGE best-fit model. The profile is reproduced to within 
3% over the entire radial range. 


TABLE 2 

MGE PARAMETERS FOR M15 


# 

log/o,v(i0Pc 

log cr'(arcmin) 

M/L 

(1) 

(2) 

(3) 

(4) 

1 

5.306 

-2.270 


2 

5.009 

-1.923 

1.9 

3 

4.822 

-1.639 

^ -^-0.5 

4 

4.670 

-1.343 

1.1 

5 

4.442 

-1.061 

1.3 

6 

4.202 

-0.784 

1 4-K).4 

7 

3.795 

-0.492 

1.7 

8 

3.361 

-0.223 

2.0 

9 

2.871 

0.022 

2.3 

10 

2.250 

0.177 

2 5 ± 0.5 

11 

2.084 

0.353 

2.5 

12 

1.127 

0.538 

2.5 

13 

1.080 

0.675 

2.5 

14 

0.066 

1.031 

2.5 


Note. — The parameters of the 15 Gaus¬ 
sians from the MGE-fit to the V-band sur¬ 
face brightness profile of Noyola & Gebhardt 
(2005). Col. (1): number of Gaussian com¬ 
ponent. Col. (2): Central surface bright¬ 
ness of each Gaussian adjusted for the as¬ 
sumed ellipticity of e = 0.05. Col. (3) Disper¬ 
sion along the major axis. Col. (4) Best-fit 
Mq/Lq^v value for each Gau ssian and as¬ 
sociated error (See Section lT^ . 

MGE fitting software developed by Cappellari (2002). Fig- 
urel^shows the comparison between the observed profile and 
the MGE fit. Table |2|gives the numerical values of the Gaus¬ 
sians that comprise the MGE model. The fit is accurate to 
better than 3% at all radii. The smallest Gaussian in the MGE 
model has a sigma of 0'^23 and the largest Gaussian has a 
sigma of 10.^7. The mass outside the tidal radius is negligible. 

We assume the photometric major axis is aligned with the 
axis of maximum rotation, as required for an axisymmetric 
model. We set the observed ellipticity of the Gaussian com¬ 
ponents to be e = 0.05, adjust the luminosity of the Gaussians 
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Fig. 5.— Observed mean velocity and velocity dispersion (first and third column) and the corresponding kinematics from the best-fit dynamical model with 
D = 10.3 kpc (second and fourth column). Top and middle row show the average kinematics of the proper motions in the x' and y'-direction and the bottom row 
shows the averaged kinematics of the line-of-sight velocities. 


accordingly to conserve flux, and take the photometric major 
axis at PA=198° north through east (Gebhardt et al. 1997). 
The chosen ellipticity sets the minimum possible value of the 
inclination of M15 to be 25°, as otherwise the Gaussians can¬ 
not be projected to have the observed flattening. VOS showed 
that modest radial variations (or errors) in the ellipticity are 
not critical for the resulting best-fit models. 

4.2. Aperture binning 

We fit the dynamical models to the observed kinematic data, 
binned in apertures on the sky. Since the models are axisym- 
metric, we first reflect all the measurements to one quadrant, 
as described in VOS, and then construct a grid of polar aper¬ 
tures. The apertures contain 50 stars per bin on average, ex¬ 
cept for the center where the bins contain only 10 stars. This 
allows an accurate measurement of the mean velocity and 
velocity dispersion (VOS). We use different sets of apertures 
for the proper motions and the line-of-sight velocities. The 
28 apertures and resulting (v) and a for the line-of-sight data 
are shown in Figure^ and given in Table 3. The line-of-sight 
apertures cover a large radial extent, from 1" to 10'. The av¬ 
erage velocity (v) and a measurement error is 1.4 km s“' and 
1.0 km s“*, respectively. The proper motion data is distributed 
over 13 apertures, also shown in Figure |3 with the corre¬ 
sponding values listed in Table 4. The average velocity (v) 
and cr measurement error is 0.04 mas yr“^ and 0.03 mas yr“^ 
respectively. The apertures extend to 18 "7 from the center. 


4.3. The parameters 

MIS is a dense and old globular cluster in which substantial 
mass segregation has taken place (e.g.. Dull et al. 1997) so 
that M/L is expected to vary with radius. Our Schwarzschild 
models therefore have not only the distance D and inclination 
i of MIS as free parameters, but must allow for a radial M/L 
variation. 

In a constant M/L model, the gravitational potential is ob¬ 
tained by multiplying the luminosity of all the Gaussians in 
the MGE mass model (Table 2) with the same M/L value. 
To construct a mass model with a smooth radial M/L profile 
we varied the M/L of the individual Gaussians as this allows 
efficient calculation of the corresponding gravitational poten¬ 
tial. However, to reduce the number of free parameters and 
to enforce a continuous profile we varied the first, third, sixth 
and tenth Gaussian and interpolated the other Gaussians log¬ 
arithmically. The Gaussians ten through fourteen were given 
the same M/L value as Gaussian ten, because their individ¬ 
ual M/L’s are not constrained well because only line-of-sight 
velocities are available at these radii. Finally, we include a 
central dark mass Mdait, represented by a point-mass poten¬ 
tial. As a result, we have seven parameters: D, i, four for the 
M/L profile, and Mdark, for which we construct models. 

4.4. Best-fit model 

Our dynamical models each have 2058 orbits covering a 
radial range of 0f'16 to 1923". Each model takes approxi- 
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Fig. 6.— Marginalized contour map of the models (crosses) varying 
inclination (horizontal axis) and distance (vertical). The three inner contours 
are drawn at formal 68.3%, 95.4% and 99.7% confidence levels for one de¬ 
grees or freedom. The subsequent contour corresponds to a factor two in¬ 
crease in Ax^. The contours outside the models are extrapolated. The best- 
fit model at inclination 60° and 10.3 kpc is denoted with a black star. The 
inclination is not well constrained. 

mately 40 minutes on an average 1.5Ghz desktop computer 
to complete. So a straightforward search of the parameters is 
impractical as this would require a minimum of 7^ = 16807 
models. Therefore we first searched the three parameters D, 
i, Mdark to find their best-fit values. After that we searched 
Mdark and the four M/L parameters keeping D, i, fixed at their 
best-fit values. Finally we checked that we found the global 
minimum by doing a small search (7^ models) through all the 
parameters. 

The model that best fits the photometric and kinematic 
observations of M15 has a total of 275 constraints and a 
= 88. To determine the error on the parameters we will 
use Ax^ = X^-xii„, where xLn is the x^ of the best fitting 
model. We will use Ax^ = 3.53 for the 68.3% confidence for 
one free parameter for D and i and Ax^ = 5.87 for the 68.3% 
confidence for 5 free parameters for Mdark and the M/L pro¬ 
file. Since our reduced x^ is much smaller than one, our use 
of Ax^ is conservative since we likely over-estimate our un¬ 
certainties. 

The best-fit model has the following parameters D= 10.3 ± 
0.4 kpc , / = 60 ± 15, Mdark = 5OO;!;5oq°M 0 and a radially vary¬ 
ing M/L profile. The best-fit M/L values are tabulated in Ta¬ 
ble 2 and shown in Figures 0 and 0 with their formal error 
bars. The best-fit kinematics are shown in Figure 0 Due to 
the small number of stars in each aperture bin the scatter in the 
observed kinematics is large. The proper motion mean veloc¬ 
ities are dominated by the errors. As a result, the mean proper 
motions of the best-fit model are very small. The overall ro¬ 
tation present in the line-of-sight velocities is fitted with the 
model. Also the observed dispersions are reproduced well. 

4.5. Inclination and Distance 

The D tan i fit of Section 13.31 gives an inclination of 
59° ±12° at a distance of 10 kpc. The contours in Figure 0 
show that the best-fitting inclination is 60° ± 15° (68.3% con¬ 
fidence, one parameters), and as such the inclination is consis¬ 
tent with the D tan i prediction, although it is not constrained 
very well. 



Fig. 7.— Radial M/L Profiles. Black line: Our best-fit V-band MGE M/L 
values and error bars signifying 68% confidence for five degrees of freedom. 
Red Line Our best-fit deprojected M/L profile. Blue solid line: The depro- 
jected M/L profile from Baumg ardt (pr iv. co mm.) Blue dashed line: Profile 
from Pasquali et al. See Section|^^and|^]for a discussion. 

The dynamical distance estimate of 10.0 ±0.5 kpc from 
M04 which assumes Ml5 is an isotropic sphere is very similar 
to our best-fit distance. Our (consistent) distance is slightly 
larger, because we allow our models the freedom to be flat¬ 
tened, and have not restricted the distribution function to be 
isotropic. 

Our best-fit distance of 10.3 ±0.4 kpc (68.3% confidence, 
one parameter) is in agreement with all other distance esti¬ 
mates; 10.4 ±0.8 kpc by Durrel & Harris (1993), 10.3 kpc 
from the globular cluster catalog of Harris (1996), 9.5 ±0.6 
kpc by Silberman & Smith (1995), and the Fe II metallicity- 
scale distance 11.2 kpc by Kraft & Ivans (2003). 

4.6. TheM/Lprofile 

Our best-fit M/L model profile is shown in black in Fig¬ 
ure 0 with the formal error bars. The error bars in Figure 0 
denote 68% confidence errors of the five parameter fit. Inter¬ 
preting the error bars is difficult, since the error bars are a one¬ 
dimensional view of the five-dimensional parameter space. 
The M/L values (and their errors) at the different radii are 
strongly correlated, since they represent enclosed mass. The 
values are listed in Table 2, and can be used to convert the 
individual Gaussian luminosity profiles into density profiles. 
Summing these provides the density profile for M15. Division 
by the luminosity profile then gives a smooth M/L profile. 
This is shown in red. 

Pasquali et al. (2004) And an Mq/Lq = 2.1 at 7' and an 
Mq/Lq = 3.7 in the center using a luminosity and mass func¬ 
tion derived from NICMOS data. Their results agree with our 
M/L profile as shown in Figure0with a blue dashed line. An¬ 
other independent measurement of the central M/L by Phin- 
ney (1993) using the acceleration of pulsars in M15 yield a 
lower limit of the Mq/Lq of 2.1 inside 0(1 and is consistent 
with our profile. 

The blue curve in Figure 0 is the radial M/L profile of 
M15 from an N-body model constructed by Baumgardt (priv. 
comm.), which is rather similar to that of Dull et al. (2003). 
It has a central peak, caused by compact remnants, mostly 
massive white dwarfs. The total number of these remnants 
that survive the cluster evolution is difficult to estimate, as 
it depends on the fraction of neutron stars that is retained in 
the cluster potential after supernova explosion, which is be¬ 
lieved to endow them with substantial kick velocities (Hansen 
& Phinney 1997; Pfahl, Rappaport & Podsiadlowski 2002). 

The profile from Baumgardt is significantly different from 
our best-fit model inside 0.'6. To be able to accurately de- 


















Dynamics of Ml5 


7 



Fig. 8. — Marginalized contour map of the models (crosses) varying the 
central Mq/Lq value (horizontal) and the dark central mass (Mq, vertical 
axis). The three inner contours are drawn at formal 68.3%, 95.4% and 99.7% 
confidence levels for three degrees of freedom. Subsequent contours con'e- 
spond to a factor 2 increase in Ax^. The best-fit model with a dark central 
mass of 500 Mq and a central M/L of 5 Mq/Lq is denoted with a black star. 
The correlation between the two parameters can clearly be seen, showing that 
the data primarily constrains the total central mass enclosed. 

termine the difference we made models with the Baumgardt 
M/L profile instead of our own. These models did not include 
a dark central mass. The Ax^ value we find for the best- 
fit model with Baumgardts M/L profile is outside our formal 
99.9% confidence level of our best-fit model. The distance 
and inclination found by using the Baumgardt profile do not 
differ significantly from our best-fit values. 

For completeness, we tried models with a constant M/L. 
The best-fit distance, inclination and central dark mass do 
not change significantly. The best-fit constant M/L found is 
1.6 ±0.2 Mq/Lq. This is the same as found by Gerssen et 
al. (2002) and consistent with 1.7 Mq /Lq found by Gebhardt 
(1997). The Ax^ value found indicates that our constant M/L 
value is consistent with our M/L profile within 95.4% (two 
sigma) confidence levels. 

Finally, the total mass of our best-fit model is 4.4 x IO^Mq. 
This is in agreement with 4.9 x 10^ Mq (Dull et al. 1997), 
4.4 X IO^Mq (M04) and 4.6 x 10^ Mq (GOO). Our estimate of 
the total mass is sensitive to the largely unconstrained value 
of M/L outside of 2', as this region contains ~ 40% of the 
mass of the cluster. 

4.7. Central dark mass 

Previous studies by Peterson et al. (1989), Gebhardt et al. 
(1997) and Gerssen et al. (2002) have argued for an interme¬ 
diate mass black hole (IMBFI) of up to a few thousand Mq 
in M15. Flowever M03 and Baumgardt et al. (2003) reported 
that they can construct N-body models for M15 which do not 
require an IMBH. For a detailed review of the history of this 
controversial subject, see M03 and van der Marel (2004). 

Our formal best-fit value of the dark central mass is 
500 / 50 °° Mq (Ax^ = 1). This mass estimate agrees with all 
the earlier estimates of the IMBH from Gebhardt (1997) to 
Gerssen et al. (2002). Figurel^shows that there is a degener¬ 
acy between the central M/L and dark central mass. 

To be able to determine the inner structure, accurate infor¬ 
mation on the central region is critical as it is difficult to de- 



Fig. 9.— Images of the central region of M15, 9" on a side. The left image 
is a WFPC2/F336W exposure. The right image is heavily massaged; we have 
convolved the original image with a boxcar of 20 pixels in size. The blue 
curves are smoothed isophotes at l75, 2'.'2 and 278 major axis radius. The 
straight line is the best-fit position angle for the isophotal major axis. 

termine what happens in the central part that is not covered by 
luminosity and kinematic data. The sphere of influence for a 
1000 Mq black hole is 0."5. The inner luminosity data point 
that we used in this model is at 0."35 from the center and the 
closest star is 0."4 and 0/25 from the center for the line-of- 
sight and proper motion data set, respectively. This results in 
a degeneracy between the central M/L and dark central mass 
as is shown in Figure |3 As a result we are not able to de¬ 
termine the nature of the matter inside ~lf'. By combining 
the mass of our best-fit mass model and the dark central mass 
we estimate that the total mass inside 1 /O arcsec (0.05 parsec, 
10“* AU at 10.3 kpc) is 3400 Mq. This implies an extremely 
high central density of 7.4 x 10“ Mq pc“^. 

To test the robustness of this density estimate we made 
models with a constant M/L (section lAhl and a dark central 
mass. In this case we found the best-fit dark central mass to be 
IOOO/Iqqq Mq. The mass contained inside 1/0 in this model 
is comparable to that from our models with the best-fit M/L 
profile. 

Guhathakurta et al. (1996) studied the cluster photometry of 
M15 using star counts. They detected 205 stars brighter than 
20th magnitude in V-band inside 1/0. The typical mass of 
these post-main-sequence and turnoff stars is 0.75 Mq, thus 
these stars account for 150 Mq. Their data is not able to con¬ 
strain the mass of fainter stars. But they give a rough (over) 
estimate and find that the total stellar mass is 4000 Mq . When 
projected onto the sky our mass model gives a total mass of 
5000 Mq in this region. This would mean that there is room 
for dark mass in the form of an IMBH with a mass in de or¬ 
der of 1000 Mq. If mass segregration has indeed taken place, 
the stellar mass function used here might significantly over¬ 
estimate the amount of low mass stars, and therefor lower the 
total stellar mass. This would allow for more dark matter in¬ 
side 1 " 0 . 

5. A DECOUPLED CORE? 

The maps of the velocity fields shown in Figure 0 display 
significant structure in the inner few arcseconds. If M15 is 
rotating, the plot of velocity against position angle for an an¬ 
nulus will show a sinuisoidal variation (it is exactly sinusoidal 
for rotation on cylinders). Figure cni shows the individual 
measurements in an annulus with inner and outer radii of Of 03 
and 0.'06 as a function of position angle on the plane of the sky 
(cf. Figure 14 of GOO). There is little azimuthal variation in 
Vx', but both Vy' and Vy' show a sinusoidal variation. Using re¬ 
lation Q, the line-of-sight velocity and y' proper motion sug¬ 
gest an inclination of 45 ±20°, consistent with our dynamical 
modeling result. Since the proper motions and radial velocites 
come from different dataset, the concordance of both results 
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Fig. 10.— Velocity of stars versus position angle on the sky for each 
velocity component. The top is for the motion parallel to the major axis, the 
middle is along the minor axis, and the bottom is the line-of-sight motion. 
The line is the best fit profile with a three-parameter fit. We fit the position 
angle of the major axis, amplitude of the rotation along minor and line-of- 
sight axes, and amplitude along the major axis. The position angle of the 
maximum in the line-of-sight velocity corresponds to the isophotal major axis 
of the central structure in M15. 

and individual significance of the rotation in each component 
strongly suggest that the central rotation in M15 is real. A 
simultaneous fit to the two kinematic data sets, subject to the 
constraint 0, sets the y'-amplitude to 11 ± 1.5 km s ' and 
jt'-amplitude to 7 ± 1.5 km s“', at a PA of 270 ± 10°. 

Figurel^shows the central region of M15 as observed with 
WFPC2/F336W. We have heavily smoothed the image on the 
right-hand side in order to measure any potential flattening. 
The image does show substantial flattening, which appears in 
a F555W image as well. The position angle of the major axis 
is not at either the kinematic PA nor at the isophotal PA of 
the major axis at large radius. We find a position angle of 
~ 120° (North to East) for the light inside of 3". This PA is 
30 degrees different from the kinematical PA as defined by the 
rotation seen in the core. Although this difference is signifi¬ 
cant given the uncertainty on the kinematically defined PA, 
the photometrically determined PA is subject to collections 
of bright stars and it is difficult to determine uncertainties for 
it. Thus, we are unable to determine whether the difference 
in the kinematic and photometric PA in the central region is 
significant. 

It is tempting to identify this misaligned structure with the 
signature of an inspiraling binary IMBFl, as described by 


Mapelli et al. (2005). Their simulations of binary IMBH pre¬ 
dict a small number of stars with high rotational velocities 
close to the black hole, and then a larger number of stars with 
aligned angular momentum (i.e., central rotation). However, 
we see no evidence for a small number of associated high ve¬ 
locity stars, which would be a direct signature in their sce¬ 
nario. We do though see ordered rotation. This rotation is 
not direct evidence for a binary black hole, but it does sug¬ 
gest an unexpected dynamical state for the central regions in 
M15. Since the relaxation time is very short in the center of 
Ml5 (around 10^ years), rotation will be quickly removed by 
two-body interaction (Akiyama & Sugimoto 1989; Kim, Lee 
& Spurzem 2004). As discussed in Gebhardt et al. (2000), it 
is difficult to maintain such strong rotation in the central parts 
of M15. A binary IMBH does increase the central relaxation 
time (by lowering the stellar density) and offers a possible 
formation mechanism—and even explains the misalignment 
between the core and halo PA—but one would need stronger 
evidence in order to invoke such a scenario. 

The misaligned core is inconsistent with our assumption of 
an axisymmetry dynamical model for M15, with major axis 
at PA of 198°. When we ignore the kinematic measurements 
inside 4" in the fitting procedure, we obtain the same distance 
and inclination. 

6. DISCUSSION AND CONCLUSIONS 

We studied the globular cluster M15 by fitting line-of-sight 
velocities, HST proper motions and surface brightness pro¬ 
files with orbit-based axisymmetric dynamical models. 

The observations used for the modeling consisted of a lumi¬ 
nosity profile from Noyola & Gebhardt (2005), 1264 line-of- 
sight velocity measurements from GOO and a sample of 703 
HST proper motions from M03. The line-of-sight data ex¬ 
tends out to 7', while the proper motions cover the inner 0.'25. 
The models provide a good fit to the observations and allow 
us to measure the distance and inclination of the cluster, the 
orbital structure, and the mass-to-light ratio M/L as a func¬ 
tion of radius. We obtain a best-fit value for the inclination of 
i = 60° ±15 and a dynamical distance of D = 10.3 ± 0.4 kpc, 
in good agreement with the canonical value. 

Our best-fit model has a 500;!;5 qq“ Mq dark central mass and 
the M/L profile shown in Figure^ which has a central peak 
and a minimum at 0.' 1. The overall shape of the profile resem¬ 
bles the shape expected for an expanded core globular cluster 
(Dull et al. 2003). The central M/L peak and the dark central 
mass together represent a mass of 3400 Mq inside the inner 
1"0 (0.05 parsec at a distance of 10.3 kpc). This suggests 
that the center harbors a large amount of dark mass. We can¬ 
not distinguish the nature of the central mass concentration. 
It could be an IMBH or it could be large number of compact 
objects, or it could be a combination. 

We found that a heavily smoothed image of M15 shows 
a flattened structure inside 4", with a different PA from the 
outer part of the cluster. The line-of-sight and proper motion 
data inside this radius can be fitted using relation (1). The 
fit gives a PA similar to that of the flattened structure and a 
rotational velocity of 10 kms“'. This suggests the structure 
is real, and constitutes a fast-spinning decoupled core at the 
center of the cluster. 

A significant improvement in the accuracy of the 
dynamical models for M15 is possible by increasing 
the accuracy of the proper motions and radial veloci¬ 
ties. This appears possible with the ACS onboard HST, 
and with high-resolution spectrographs on 8m class tele- 
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TABLE 3 

Kinematics of the proper motions in polar apertures 


# 

(1) 

n* 

(2) 

^0 

(3) 

do 

(4) 

Ar 

(5) 

A9 

(6) 

U' 

(7) 

AV,, 

(8) 

cr,, 

(9) 

Act,, 

(10) 

V,, 

(11) 

AV,, 

(12) 

o-,, 

(13) 

Act,, 

(14) 

1 

50 

0.021 

0.785 

0.033 

1.571 

-0.04 

0.04 

0.25 

0.03 

-0.02 

0.04 

0.24 

0.04 

2 

51 

0.060 

0.393 

0.046 

0.785 

-0.02 

0.04 

0.21 

0.02 

-0.02 

0.04 

0.22 

0.03 

3 

54 

0.060 

1.178 

0.046 

0.785 

-0.04 

0.03 

0.23 

0.02 

-0.01 

0.04 

0.25 

0.03 

4 

54 

0.108 

0.262 

0.050 

0.524 

0.01 

0.04 

0.28 

0.02 

-0.00 

0.03 

0.19 

0.03 

5 

54 

0.108 

0.785 

0.050 

0.524 

0.06 

0.04 

0.22 

0.03 

0.07 

0.03 

0.19 

0.03 

6 

63 

0.108 

1.309 

0.050 

0.524 

0.02 

0.03 

0.23 

0.03 

0.04 

0.03 

0.17 

0.03 

7 

45 

0.154 

0.262 

0.043 

0.524 

0.02 

0.04 

0.22 

0.02 

0.04 

0.04 

0.26 

0.03 

8 

48 

0.154 

0.785 

0.043 

0.524 

-0.06 

0.03 

0.21 

0.03 

-0.01 

0.03 

0.15 

0.03 

9 

61 

0.154 

1.309 

0.043 

0.524 

0.03 

0.03 

0.19 

0.02 

-0.02 

0.02 

0.19 

0.02 

10 

56 

0.214 

0.262 

0.077 

0.524 

-0.02 

0.02 

0.15 

0.03 

-0.02 

0.03 

0.22 

0.02 

11 

58 

0.214 

0.785 

0.077 

0.524 

0.02 

0.03 

0.21 

0.03 

0.00 

0.03 

0.25 

0.03 

12 

62 

0.214 

1.309 

0.077 

0.524 

-0.05 

0.03 

0.18 

0.02 

-0.01 

0.03 

0.25 

0.03 

13 

47 

0.279 

0.785 

0.053 

1.571 

0.01 

0.04 

0.21 

0.02 

-0.03 

0.04 

0.24 

0.02 


Note. — The mean velocity and velocity dispersion of the proper motion observations calculated in polar apertures on the plane of the 
sky. Per row the information per aperture is given. The first column labels the aperture and the second column gives the number of stars n* 
that fall in the aperture. Columns 3—6 list the polar coordinates r (in arcmin) and the angle 9 (in degrees) of the centroid of the aperture and 
the corresponding widths Ar (in arcmin) and A0 (in degrees). The remaining columns present the average proper motion kinematics in units 
of masyr^*. The mean velocity V with error AV and velocity dispersion a with error Act are given in columns 7—10 for the proper motion 
component in the x'-direction and in columns 11-14 for the proper motion component in the y'-direction. 
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TABLE 4 

Kinematics of the line-of-sight velocities in polar apertures 


# 

(1) 

(2) 

'•0 

(3) 

00 

(4) 

Ar 

(5) 

Ae 

(6) 

V,, 

(V) 

AV., 

(8) 

CT,/ 

(9) 

Act-/ 

(10) 

I 

10 

0.015 

0.785 

0.022 

1.571 

0.6 

3.8 

10.8 

2.6 

2 

10 

0.030 

0.785 

0.008 

1.571 

1.0 

4.1 

10.4 

2.5 

3 

10 

0.040 

0.785 

0.011 

1.571 

-0.7 

4.8 

13.0 

3.8 

4 

10 

0.049 

0.785 

0.007 

1.571 

-8.9 

4.2 

11.2 

2.2 

5 

62 

0.103 

0.393 

0.102 

0.785 

0.2 

1.6 

11.3 

1.1 

6 

56 

0.103 

1.178 

0.102 

0.785 

1.2 

1.8 

12.8 

1.4 

7 

55 

0.214 

0.262 

0.120 

0.524 

1.4 

1.4 

9.4 

0.9 

8 

63 

0.214 

0.785 

0.120 

0.524 

2.3 

1.3 

9.3 

1.0 

9 

52 

0.214 

1.309 

0.120 

0.524 

3.2 

1.4 

9.4 

0.9 

10 

63 

0.366 

0.196 

0.184 

0.393 

0.6 

1.4 

10.3 

1.0 

11 

63 

0.366 

0.589 

0.184 

0.393 

0.5 

1.4 

10.3 

1.0 

12 

54 

0.366 

0.982 

0.184 

0.393 

1.8 

1.9 

12.3 

1.1 

13 

70 

0.366 

1.374 

0.184 

0.393 

1.7 

1.2 

10.2 

0.9 

14 

53 

0.553 

0.196 

0.190 

0.393 

0.7 

1.5 

10.2 

1.2 

15 

57 

0.553 

0.589 

0.190 

0.393 

3.6 

1.5 

9.8 

0.8 

16 

63 

0.553 

0.982 

0.190 

0.393 

2.5 

1.2 

9.6 

0.8 

17 

60 

0.553 

1.374 

0.190 

0.393 

0.7 

1.3 

9.5 

0.8 

18 

61 

0.791 

0.196 

0.286 

0.393 

2.7 

1.3 

9.0 

0.8 

19 

69 

0.791 

0.589 

0.286 

0.393 

2.5 

1.2 

9.3 

0.8 

20 

78 

0.791 

0.982 

0.286 

0.393 

2.9 

1.1 

8.5 

1.0 

21 

48 

0.791 

1.374 

0.286 

0.393 

3.4 

1.5 

9.9 

2.1 

22 

68 

1.165 

0.196 

0.462 

0.393 

3.4 

1.0 

7.8 

0.8 

23 

49 

1.165 

0.589 

0.462 

0.393 

1.0 

1.3 

8.0 

1.2 

24 

66 

1.165 

0.982 

0.462 

0.393 

0.6 

1.1 

8.5 

1.4 

25 

62 

1.165 

1.374 

0.462 

0.393 

1.2 

1.2 

9.1 

2.0 

26 

60 

2.560 

0.393 

2.329 

0.785 

1.8 

0.8 

6.5 

0.7 

27 

66 

2.560 

1.178 

2.329 

0.785 

0.4 

0.9 

6.7 

1.0 

28 

60 

5.306 

0.785 

3.163 

1.571 

1.5 

0.6 

3.9 

0.4 


Note. — The mean velocity and velocity dispersion calculated in polar apertures on the plane of sky from the 
line-of-sight velocity observations. Columns 1-6 are as in Table|^and the remaining columns present the average 
line-of-sight kinematics in km s“*. 






